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Abstract. The 7.8 eV nuclear isomer transition in ^^^Thorium has been suggested 
as a clock transition in a new type of optical frequency standard. Here we discuss 
the construction of a "solid-state nuclear clock" from Thorium nuclei implanted into 
single crystals transparent in the vacuum ultraviolet range. We investigate crystal- 
induced line shifts and broadening effects for the specific system of Calcium fluoride. 
At liquid Nitrogen temperatures, the clock performance will be limited by decoherence 
due to magnetic coupling of the Thorium nuclei to neighboring nuclear moments, 
ruling out the commonly used Rabi or Ramsey interrogation schemes. We propose 
clock stabilization based on a flourescence spectroscopy method and present optimized 
operation parameters. Taking advantage of the high number of quantum oscillators 
under continuous interrogation, a fractional instability level of 10~^^ might be reached 
within the solid-state approach. 
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1. Introduction 

Atomic clocks liave a significant impact both on fundamental science as well as 
on everyday life technology. Prominent examples for the latter are satellite-based 
navigation systems, such as GPS, GLONASS, and GALILEO, and telecommunication 
networks. On the fundamental side, the current definition of the unit of time (SI 
second) is a duration corresponding to 9,192,631,770 periods of microwave (MW) 
radiation stablized to the ground state hyperfine transition in Cesium [1]. The 
physical implementation is a global network of primary Cesium frequency standards, Cs 
atomic fountain clocks have reached an uncertainty level of 4.5 ■ 10~^^ [2]. In fountain 
clocks the interrogation time is limited by the ballistic flight time of the atoms within 
the apparatus, further progress within this approach can be achieved in space-based 
realizations as the PHARAO project [3]. 

The recent progress in atom cooling and trapping together with the emergence of 
sub-hertz linewidth lasers allows to develop more accurate frequency standards based 
on optical transitions using trapped neutral atoms and ions. For example, an inaccuracy 
level evaluation for a neutral ^^Sr optical lattice clock is 1.5 ■ 10~^^ [4], and the ^^Al"*" 
quantum logic frequency standard demonstrated an inaccuracy level of the order of 
10~^^ 0. These approaches benefit from a 10^ — 10^ gain in quality factor by going from 
microwave to optical transitions and longer observation times using trapped quantum 
oscillators. 

A promising candidate for the role of a clock transition in a next-generation 
quantum frequency standard is the nuclear isomer transition in ^^^Thorium. This unique 
radio isotope has an extremely low-energy 229mrj]^ isomeric state. The most recent 
estimate of the transition energy is 7.8 ± 0.6 eV, obtained from indirect measurements 
of the ^^^U— J-^^^Th decay spectrum with an advanced X-ray microcalorimeter [6]. The 
corresponding wavelength of about 160 nm in the vacuum ultraviolet (VUV) is attainable 
by modern laser systems and hence the construction of a Thorium "nuclear clock" seems 
plausible. 

An optical frequency standard based on a nuclear transition will open new 
possibilities for experimental studies of basic laws of physics, i.e. constraining possible 
drifts of fundamental constants of interactions. In [7] it is conjectured that a possible 
drift of the fine structure constant 6a/ a will result a corresponding but amplified change 
of the ^^^Th nuclear transition frequency 600/00 with an amplification factor between 
10^ and 10^. Hence the Thorium quantum frequency standard will be several order 
of magnitude more sensitive to possible drifts of a than existing standards based on 
electronic shell transitions in atoms or ions. 

It should be noted that the current energy value of the Thorium isomer transition is 
derived indirectly by looking at 7 transitions with energies much higher than the isomer 
one. Therefore this value can not be considered definitive until a direct measurement 
is performed. Several groups worldwide (i.e. PTE, Germany; Georgia Institute of 
Technology, USA; University of California, Los Angeles, USA; Vienna University of 



Performance of a Thorium solid-state nuclear clock 



3 



Technology, Austria) are currently working towards this goal. 

Two main experimental approaches towards a Thorium clock, first proposed in [8], 
are under realization now: the use of trapped Thorium ions [9]-[TTj and a "solid-state 
nuclear spectroscopy" [I2|[T3]. 

In the first scheme, individual or small ensembles (up to 10^) of ions are trapped 
and laser cooled using established techniques. The main advantage of the trapped ion 
approach is a precise knowledge of and control over the fields affecting the ions. It has 
been proposed that a nuclear clock based on a single trapped 229'j]^3+ might reach 
a fractional inaccuracy of 10^^^ |17j . 

In the second approach, a macroscopic number (up to 10^^) of Thorium atoms 
is implanted into a VUV-transparent crystal and interrogated simultaneously. This 
scheme will require none of the complex experimental techniques for ion trapping and 
cooling and can be fabricated in mass production. The number of simultaneously excited 
Thorium nuclei is many orders of magnitudes higher than in the trapped ion approach, 
therefore the solid-state spectroscopy method seems to be more promising for an initial 
localization of frequency of the isomer transition. 

In the present work we analyze the possibility of constructing a "solid-state nuclear 
clock" based on Thorium-doped crystals. Here, certain difficulties arise from the 
systematic crystal effects and the very long relaxation time of the excited nuclear state 
that can significantly exceed the interrogation time. We briefly review the current state 
of knowledge on the Thorium isomer energy and lifetime and consider transition line 
shifts, splittings, and broadening processes due to the crystal environment for the specific 
case of Thorium doped into Calcium flouride (CaF2) crystals (section [2]). We then 
describe the interaction of the nuclear quantum states of the Thorium doping complex 
with coherent laser radiation in section [31 In section H] we propose an interrogation 
scheme based on counting of fluorescence photons to implement a solid-state nuclear 
clock and estimate the attainable fundamental limits of precision. At liquid Nitrogen 
temperatures, we find the clock performance to be determined by spin-spin relaxation 
due to interactions of the Thorium nuclear moment with neighboring nuclei. A detailed 
analysis of this process (presented in section [5]) for the Thorium-doped CaF2 system 
indicates, that a fractional instability on the 10~^^ level might be reached (neglecting 
technical limitations i.e. imposed by the interrogation laser system) within the solid- 
state clock approach. 

2. Thorium nuclear transition in the crystal lattice environment 
2.1. Transition energy and isomer state lifetime 

The isomer transition energy of ^^^Th is currently derived from indirect measurements 
of the 7-ray spectrum resulting from the decay of ^'^^Uranium. In 1989-1993 Helmer 
and Reich performed first measurements using high quality Germanium detectors 
(resolution from 300 to 900 eV). They estimated the ^^^Th isomer state energy to 
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he E = 3.5 ± 1.0 eV [TSl[T2j. This unnaturally low value triggered a multitude of 
investigations, both theoretical and experimental, trying to determine the transition 
energy precisely, and to specify other properties of the isomer state of ^^^Th (such as 
the lifetime and the magnetic moment). However, searches for direct photon emission 
from the low- lying excited state have failed to report an unambiguous signal [20H23] . In 
2005, Guimaraes-Filho and Helene re-analysed the old experimental data and reported 
the value oi E = 5.5 ± 1.0 eV [21]. New indirect measurements with an advanced X- 
ray micro calorimeter (resolution from 26 to 30 eV) were performed by Beck et al. in 
2007 [25]. They published a new value for the transition energy E = 7.6 ± 0.5 eV for 
the isomer nuclear transition, shifting it into the vacuum ultraviolet domain. This shift 
probably explains the absence of signatures of the transition in previous experiments. In 
2009, Beck et al. reviewed their results taking into account the non-zero probability of 
the 42.43 keV 229mrp]^ transition (estimated as 2 %) and published a revised version 
E = 7.8 ± 0.5 eV [6]. This value is the currently most accepted one in the community 
but can not be considered definite until a direct measurement is performed successfully. 
The dominant uncertainty in the prediction of Beck et al. is connected with the value 
of the branching ratio b = 1/13 [25] from the 29.19 keV level to the ground state [26] . 
Estimation of this value vary in different references. In [26] two alternative values for 
this branching ratio are mentioned: 25 % which would result in E = 9.3 ± 0.6 eV 
and 51% which would result in E = 14.0 ± 1.0 eV. In the present analysis we use 
the value of 7.8 ± 0.5 eV, scaling the obtained results to different transition energies is 
straightforward. 

While a direct measurement of the isomer transition energy in ^^^Th is outstanding, 
the decay rate 7 of the isomer state also remains unknow. The most recent estimations 
of the half life of the isomer state for a bare nucleus are based on theoretical calculations 
of the matrix element of the magnetic moment [27]. The theory was verified by 
comparison with experimental data for transitions at higher energies. It predicts a 
half-life of T1/2 = (10.95h)/(0.025-E'^) for the isomer transition, where E is given in eV. 
For E = 7.8 eV this yields T1/2 = 55min which corresponds to a spontaneous decay rate 
of 7s = 2.09-10"^ s~^. In [28] it is predicted that the probability of spontaneous magnetic 
dipole emission in a transparent non-magnetic dielectric medium with refractive index 
n is times higher than in vacuum. The reason for this factor is the renormalisation of 
the density of emitted photon states. For a CaF2 crystal, n = 1.55 at a wavelength 
A = 160 nm (corresponding to E = 7.8 eV), the enlargement factor is 3.7 for the 
spontaneous decay rate of the isomeric state. Therefore, if the crystal contains 
excited Thorium nuclei, the number dNph of 160 nm photons emitted during a short 
time interval dt can be estimated as 

dNph = N.^sn^dt. (1) 

The total decay rate 7 is determined not only by spontaneous decay of the nucleus, it 
can be enhanced by interactions between the nucleus and the electron shell in electronic 
bridge processes and bound internal conversion effects [29l[3Tl[32]. Such a decay can 
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be accompanied by the emission of several low-energy photons with widely distributed 
energies [29]. This would make it difficult to separate decay events from optical noise 
sources (i.e. from scintillation or radioluminescence). 

There are a number of theoretical calculations on electronic bridge processes in 
the Thorium system. In an isolated neutral atom the Thorium isomer state lifetime is 
expected to be about 4.5 min for E = 3.5 eV and about 10~^s for E = 7.6 eV [29] 0. 
In the Th^'^ ion for E = 7.6 eV the expected lifetime is the same as in a bare nucleus 
described above, if the valence electron is in the ground 5/5/2 state. It decreases by 
a factor of 20 if the valence electron is in the metastable 7s state [31]. In the T/i+ 
ion the lifetime decreases 10^ — 10^ times for E = 3.5 and E = 5.5 eV [32], there is 
not sufficient data about the Th"*" ion electronic structure to calculate the lifetime at 
E = 7.8 eV. Similar calculations for the r/i^"*" ion, as it will occur in the solid-state 
approach presented here, have yet to be performed. 

There were several experimental attempts to determine the ^"^^"^Th lifetime without 
knowledge of the exact isomer energy. Inamura et al. measured a temporal variation of 
the a decay activity in a sample of ^^^Th after excitation in a hollow-cathode electric 
discharge [33]. They propose a value of T1/2 = 2 ± 1 min for the half-life of the isomeric 
state which corresponds to a total excited state decay rate of 7 = 6^2 ' 10~^s~^. 
Kikunaga et al. studied the a spectrum of chemically "fresh" ^^^Th produced in 
the a decay of ^ss-y r^^ie authors obtain T1/2 < 2h with 3a confidence. The 

latter experiment was not suited to measure half-lives on the order of minutes because 
of the long chemical preparation process. The analysis of both these experiments is 
significantly complicated by the presence of various chemical compositions of Thorium, 
mainly ThO''', Th and Th"*" in [33] and hydroxide and chloride complexes in [34]. These 
chemical structures determine the energies of electronic levels and hence the rates of 
electronic bridge processes. Therefore the different results obtained in [331 [M] may 
be explained by the different chemical structures and the times of preparation and 
measurement. 

In the present work we assume the band gap of the Thorium-doped CaF2 material 
to be sufficiently large to neglect electronic bridge processes and contributions from 
low-energy photons to the measured fluorescence count rates. In other words, we expect 
no electronic levels near 7.8 eV in the Thorium-doped CaF2 crystal. This assumption 
will be substantiated in a future publication [16]. We perform our analysis using the 
theoretical value for the bare nucleus lifetime j = jg-n^ = 7.8-10~^s~^, enhanced by the 
refractive index, and assume relaxation via electronic bridge processes to be negligible. 
Again, as for the energy, scaling our results to different values of the 229mrj]^ lifetime is 
straightforward. 

I Thorium first, second, third, and fourth ionization energies are 6.31 eV [3D], ll.SeV, 20eV, and 29eV 
respectively. Therefore the expected energy of the isomer transition (7.8 eV) is above the ionization 
energy of the neutral Thorium atom. 
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2.2. Effects of the crystal lattice environment 

In this work we consider the specific system of Thorium nuclei doped into Calcium 
fluoride single crystals. CaF2 is a standard material in UV optics, it is routinely grown 
and doped in many laboratories and companies worldwide using various techniques 
(Czochralski, Bridgman-Stockbarger, micro-pulling-down). Standard procedures for 
cutting, cleaving, polishing, and coating are established, CaF2 is well characterized 
concerning radiation damage, both concerning UV laser radiation and radiation related 
to the nuclear decay chain of ^^^Th. 

CaF2 has a wide bandgap of about 12 eV, values in the literature range from 
11.6 eV [11] to 12.1 eV [I5]. The valence band of CaF2 consists of 2p levels of the 
Fluorine ions, while the bottom of the conduction band originates from the 4s and 3d 
orbitals of the Calcium ions. Optical transmission of better than 70 % in the 160 nm 
wavelength regime is routinely achieved in high-purity single crystals. 

The simple cubic lattice structure and the fact that only the Ca^"*" site qualifies 
as a doping site makes Thorium-doped CaF2 amendable to full electronic structure 
simulations [16] which is currently out of reach for significantly more complex crystal 
structures (i.e. LiCaAlFg as proposed in [T2]). 

Thorium atoms in a solid-state crystal lattice are confined in the Lamb-Dicke 
regime, the recoil energy Erec = (7.8 eV)^/(2Mr/iC^) ~ 1.4 ■ 10"^° eV is far below the 
energy required to create a phonon with a recoil momentum Prec = 7.8 eV/c in the lattice. 
Therefore the only process that could transfer the recoil energy to the lattice vibrations 
must include the scattering of optical phonons. However optical modes are expected to 
be frozen out even at room temperature [12]. Therefore, internal and external degrees 
of freedom are decoupled and there is no sensitivity of the nuclear transition to recoil 
or first-order Doppler effects. However, the interaction between the Thorium nuclear 
levels and the crystal lattice environment can cause some inhomogeneous effects which 
will affect the performance of a solid-state nuclear clock. 

We assume the ^^^Th ions to be chemically bound into a VUV transparent ionic 
crystal. Therefore the fine interaction of the Thorium nuclei with the crystal lattice 
environment is absent because of the absence of unpaired electronic spins. The remaining 
interaction is the hyperfine interaction where the Hamiltonian H' = Hrfs can be 
represented in a multipole expansion Hhfs = Heq + Hmi + He2 + ••• [III- In the 
following we will discuss the individual terms in more detail. 

The electric monopole shift {Heo) arises from contact interaction between the 
electron cloud and the finite nuclear volume. As the ground state and the excited nuclear 
state can have different sizes, this effect leads to a shift of the nuclear transition relative 
to a bare nucleus. Typical monopole shifts are about 1 GHz [35], they will depend 
on the specific choice of the host crystal. However, this effect will be identical (up to 
temperature effects) for all nuclei having the same position within the lattice and hence 
deteriorate the solid-state nuclear clock accuracy, but not the stability. Temperature 
effects appear because the electron density at a given lattice site is temperature 
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Figure 1. CaF2 crystal lattice (a) and Thorium inclusion with charge compensation 
by two in-line Fluorine interstitials (b). 

dependent. This temperature dependence is estimated as lOkHz/K [12] which leads 
to a broadening of about 10 Hz for a temperature stabihty of AT ^ 1 mK which is 
attainable by standard techniques. It should be noted that a certain contribution to 
the decoherence rate can be caused by crystal lattice vibrations. These vibrations will 
modulate the electronic density in the vicinity of the nucleus thus broaden the transition 
line. A detailed study of this broadening effect is beyond the scope of this article. We 
expect that this effect, if relevant at all, can be effectively suppressed by cooling of the 
crystal lattice. 

The electric quadrupole shift {He2) appears due to the electric quadrupole moment 
Q of the ^^^Th nucleus in the presence of a non-vanishing electric field gradient at 
the position of the nucleus. Reported experimental values for the quadrupole moment 
in the ground state are Qg = 3.149(32) eb [36] and Qg = 3.11(16) eb [37], the value 
Qis of the quadrupole moment in the isomer state is still unknown, estimations [38] 
give Qis ~ 1.8 eb. If the electric field gradient in the principal axes has two equal 
components, 4)xx = <Pyy = —'Pzz/'2, the matrix element of the quadrupole interaction 
term can be written as [39] 

{Im\HE2\Im') = Smm'^j^^^ [3m' - I{I + 1)] (2) 

where the z axis is the axis of symmetry of the electric field, (pzz = d'^^p/dz', I and m are 
the nuclear angular moment and its projection on the axis of symmetry of the electric 
field respectively. In a Th-doped CaF2 crystal, the Thorium ion will most probably 
replace a Ca^"*" ion in a Th^"*" state. To recover an ionic (and hence transparent) crystal, 
the two remaining electrons will be localized by a charge compensation mechanism, 
either by introducing two additional Fluorine (F~) interstitials or by creating a vacancy 
by leaving a neighboring Ca^"*" site empty. First simulations indicate a preference for 
charge compensation by interstitials which will be the main contribution to the electric 
field gradient [IQ]. 

Consider the charge compensation by two F~ interstitials as shown in Figure [T](b). 
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Figure 2. Quadrupole structure of the ground and isomer state of the Th nucleus 
in tlie CaF2 crystal lattice with charge compensation as presented in Figure [Ijb). 
Light arrows represent the possible magnetic dipole transitions, matrix elements are 
proportional to the indicated Clebsch-Gordan coefficients C/^"^ = m' — m). The 
quantization axis z is chosen along the axis of symmetry of the electric field gradient. 



We expect the geometrical lattice distortion around the dopant to be negligible. 
Using the Ca^"*" lattice constant a = 5.6 A and the quadrupole antishielding factor 
7oo = —177.5 IH] we estimate ~ —5.1 ■ lO^^V/cm^. This leads to quadrupole 
splittings of the order of a few hundred MHz which can be easily resolved in laser 
spectroscopy (see Figure [2]). Note that this splitting could also be directly measurend 
experimentally using standard NMR techniques, which would give valuable information 
on the details of the charge compensation mechanism. 

The quadrupole interaction of the nucleus with phononic modes of the crystal lattice 
leads to inhomogeneous broadening and mixing between different quadrupole sublevels 
as the the typical distances between the quadrupol sublevels are much smaller than 
characteristic phonon energies. Comparison with experimental data for the quadrupole 
spin-lattice relaxation of the Halogen nuclei in CsBr and Csl crystals |12] and scaling 
the theory of Van Kranendonk |13] and Yoshida and Moriya [S] to the CaF2 crystal 
yields a relaxation and mixing rate of the order of 27r ■ (5 — 100) Hz at room temperature. 
This rate scales as T^, leading to about 2tt ■ (0.3 — 7) Hz at liquid Nitrogen temperature. 
The temperature of the crystal should however not be so low to entirely suppress the 
mixing, as populations accumulating in quadrupole sublevels not interacting with the 
laser field could lead to a loss of signal or neccessitate additional repump lasers. 

The 2^'^ order Doppler effect leads to a shift and broadening of the isomer transition 
due to vibrations of the nuclei. Using a classical harmonic oscillator model, we 
find the shift S = -u{v'^) /{2c^) = -3ukBT/{2MThc'^) and the broadening Au = 
oj (Av^) / {2c'^) = ^/QuksT / {2MThC^)- Here oj is the isomer transition frequency, v the 
velocity of the Thorium nuclei, (Af ^) = ^J {v^) — {v"^)"^. This leads to a shift of about 
271 ■ 1.14 • T (in Hz), and to a broadening of about 27r ■ 0.93 ■ T (in Hz) where T is 
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the temperature in Kelvin. At liquid Nitrogen temperature we expect a shift of about 
2tt ■ 88 Hz and a broadening of about 2tt ■ 72 Hz. 

The term Hmi, i-e. the magnetic dipole interaction with randomly oriented spins 
of surrounding nuclei also contributes to broadening. The characteristic energies of 
spin-spin interactions correspond to the temperature range of the order 1 microKelvin, 
the broadening is expected to be temperature-independent for experimentally relevant 
conditions, in contrast to all other broadenings discussed above, which decrease with 
temperature. A detailed evaluation (see Appendix A) yields a contribution to the optical 
relaxation rate from ~ 27r ■ 84Hz for the |/ = 5/2, m = ±3/2) ^ |/' = 3/2, m' = ±1/2) 
optical transition to ~ 27r-251Hz for the |/ = 5/2, m = ±5/2) ^ |/' = 3/2, m' = ±3/2) 
optical transition. 

To resume, we expect a global uniform shift of about 1 GHz dominated by the 
electric monopole shift compared to the bare nucleaus isomer transition frequency 
(assuming identical doping complexes and charge compensation). Main sources of 
inhomogeneous broadening are magnetic dipole interaction and 2"*^ order Doppler effect, 
leading to an ensemble linewidth of 2tx ■ (84 — 251) Hz and 27r ■ 0.93 Hz ■ T respectively 
where T is the temperature in Kelvin. Taking the ensemble linewidth of 27r ■ 150 Hz we 
expect a quality factor Q = 1.26 ■ 10^'^ for the solid state clock transition. Note that 
crystal field shifts and broadenings (and their temperature dependence) can directly be 
studied with narrow band lasers, transferring Mofibauer spectroscopy into the optical 
regime. 

3. Excitation of Thorium nuclei by laser radiation 

Let us consider an ensemble of 229rji]-^4+ ^^^^ implanted into a CaF2 crystal and 
interacting with a narrow-band vacuum ultraviolet laser whose frequency Ul is close 
to the frequency u of the isomer transition. Here we consider only the ions in one 
charge-compensating configuration, i.e. the one presented in Figure [T](b). All other 
configurations will not interact with the laser field due to different quadrupole shifts. 
Th^"*" ions are iso-electronic with the noble gas Radon with zero angular momentum 
and zero spin due to the electronic shell. In the following we neglect the electronic shell 
structure as well as higher excitation levels of the Th nucleus. The nuclear ground and 
isomer states are degenerate with respect to the sign of projections m and m' and with 
respect to angular momenta / and I' up to interactions with randomly oriented spins 
of surrounding nuclei. The magnetic field B of the spectroscopy laser radiation can be 
written a.s B = Bq cosw^t. 

The Hamiltonian H of the Thorium nucleus in resonance approximation can be 
represented as 

H = Ho + hY,Veg 0^7)(e|e'""* ± |e) ((^le"*-^*] ± H' . (3) 

Here, the indices e and g indicate the sublevels of the nuclear excited and ground 
state respectively, V^g = B^qulg/ {2h), BQg is the q-th covariant component of Bq, /if^ 
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is the q-th contravariant component of the matrix element /leg of the nucleus magnetic 
momentum operator in cyclic coordinates, q = m' — m, Hq is the Hamiltonian of the 
free Thorium ion. The term H' describes the interactions of the nucleus with the crystal 
lattice leading to shifts (which could also be added to Hq) and relaxation terms which 
will be discussed later. The set of equations for the density matrix pij with relaxation 
terms in the rotating frame reads: 

Pee = - ^ Kg [pge - Peg] - IPee + Pe'e'] (4) 

a e' 

Pee' = - ^ l^egPge' " K'gPeg] " (^^ee' + ^±)Pee', 6 c' ] (5) 

9 

^i{uJL - UJeg) " T' Peg - ^eg'Pg'g + i ^ Ve'gPee'] (6) 



Peg 
P99 



-i^Veg [peg " Pge] + legpee + ^ Pg'g' ; (7) 

e eg' 

P99' = -^Y^ i^egPeg' - V^g'Pge] " {iujgg' + T^^)pgg', Q 7^ q' ■ (8) 
e 

Here F^^ and T^^^^ are decoherence and mixing rates between magnetic states of the 
same quadrupole energy level, V is the relaxation rate of optical coherence peg, ujij is 
the difference between the frequencies of the i-th and j-th level, 7eg = 7 ■ |C*/^'gp is the 
rate of spontaneous |e) — > \g) transitions. 

In this work we suppose that the spectroscopy laser field is weak, i.e. 

Kg«min[r',r^^r|]'^',rf']. (9) 

Then, from equations (H]) - (j8]) and from the normalization condition pa = 1 we 
derive: 

^gg'iX ^ Pexc) r pexc „n 

P99' ~ 2/ + 1 ' ~ 21' + 1 ' ^ ' 

where pexc = XlePee IS the total excited state population. From equations ([6]), ([8]), and 
ffTU]) we obtain: 

E2r'Vgg /I hpexc\ 

r'2 + (a;^-a;,J2 \l~ ^ ) " (^^) 

The quadrupole splitting in the ground and the isomer state is of the order of hundred 
MHz which is much greater than the expected value of V (about a few kHz). Hence 
only transitions between certain quadrupole sublevels of the ground and excited states 
should be taken into account. Let us consider the case where then laser is tuned close 
to the transition between \g1.2) = 1-^ = 5/2, m = ±3/2) and lei.2) = \I' = 3/2, m = 
±1/2). This is the transition with the lowest relaxation rate of optical coherence, see 
Appendix A. In the following we will name this transition the clock transition and use 
Ueg = oj. We introduce the optical detuning /S. = ul — uj. Suppose, for the sake of 
simplicity, that the radiation has a certain polarization along g'-th cyclic unit vector, so 
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^S(Ml)(2/' + 1) = 

0.65 /^Tv is a reduced dipole moment (i?(Ml) = 0.025 /i^ is the reduced transition 
probabihty [27]), Cf'^^ is a Clebsch-Gordan coefficient, c is the speed of hght in vacuum. 
Using J2m,q |C/Xr = 1) we express 
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Averaging over possible spatial orientations of the electric field gradient with respect to 
the polarization of the laser, and expressing the amplitude \Bq\ via the intensity Jq of 
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Here we denote the averaging over spatial orientations by angled brackets and we have 
taken into account the values of the ^^^Thorium nuclear angular momenta 7 = 5/2 and 
/' = 3/2 in the ground and excited state respectively. 



4. Interrogation scheme and fundamental limits of stability 

Any passive quantum frequency standard consists on three main elements: a generator 
producing some periodic signal (a laser in the case of an optical frequency standard) 
with good short-time stability and a frequency that can be tuned around a "nominal" 
frequency Un, a quantum discriminator, and a servo loop determining and correcting the 
offset f = LOn — OJ between the nominal frequency cj„ and the frequency u) of the clock 
transition according to the result of the interrogation of the quantum discriminator. If 
we add a system to count the number of oscillations of the generator, we obtain a clock 
(see Figure [3]). 

Usually, the Rabi, Ramsey or hyper- Ramsey [IS] interrogation scheme is used 
in high-performance frequency standards. However, all these methods require the exact 
knowledge of the quantum state of the quantum discriminator at the beginning of any 
interrogation cycle. If the quantum discriminator consists of one or a small number of 
ions it is possible to know the quantum state precisely at the end of any interrogation 
cycle (and therefore at the beginning of any next cycle) as the interrogation measurement 
leads to a projection of the quantum state onto a well-defined basis. If the quantum 
discriminator consists of a large number of trapped ions it is possible either to reload 
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discriminator 




Figure 3. Schematic representation of a passive quantum frequency standard 
constructed from a generator (laser), a quantum discriminator, and a feedback loop. 
Linking the stabilized generator to an electronic counting system (i.e. through a 
frequency comb system) allows the realization of a clock. 

the trap with ions of a known quantum state or to repump the ions back to the initial 
quantum state via some short-hving state. Note that for Thorium ions in an ion trap 
such "repumping" of the isomer state could be performed using an excitation of the 
electronic shell to a resonant metastable state which will lead to an amplification of 
electronic bridge processes [3Tj . 

In the solid-state nuclear frequency standard such a repumping seems to be 
impossible due to the absence of an accessible metastable state in the band gap. Another 
fundamental problem arised from the huge difference between the relaxation time of the 
excited state population (presumably tens of minutes) and of the short coherence time 
between the ground and excited state (milliseconds) due to crystal lattice effects. At 
least one Rabi or Ramsey interrogation cycle should be completed during the time 
interval t < l/V whereas the time tp needed for bringing the nuclei back to the ground 
state should be > I/7. Therefore the Rabi or Ramsey interrogation time would take 
only a tiny (f» 10~^) fraction of the total interrogation cycle which would significantly 
reduce the clock performance due to the Dick effect [17] . Moreover, Ramsey or Rabi 
interrogation schemes require quite strong laser fields, able to excite a significant number 
of nuclei during a short interrogation time. 

4.I. Clock interrogation by fluorescence spectroscopy 

To tackle these problems we propose an interrogation scheme based on counting of 
spontaneous nuclear decay fluorescence photons after illumination of the quantum 
discriminator by the spectroscopy laser. The (n, n + 1) interrogation cycle consists 
of four time intervals, as depicted in Figure IH In the first interval 6 between tn and t'„ 
the crystal sample is illuminated by laser radiation with the frequency detuned by 6m 
to the blue side of the nominal frequency Un of the generator. Then a laser shutter is 
closed and the fluorescence photons are counted by a photodetector during the second 
time interval 9' between t'^ and tn+i- Then these two operation are repeated, but with 
the laser frequency shifted to the red side by 6m of the nominal value. After the last 
measurement we determine the frequency offset f = Un — u and make a correction to the 
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Figure 4. Interrogation cycle of quantum discriminator as discussed in the text. Note 
that this cycle is opperatcd continuously without any deadtime. 



nominal frequency. As no further state initialization, repumping, or sample preparation 
is required, the next interrogation cycle can start imediately. 

We will now consider a single interrogation cycle in more detail. Solving 
equation we express the excited state population p^xc at time t'^ as 

5R/2 



Pexcifn) — Pexcitn) GXp 
+ 



R 



i + {5m + f)yT'- 



l[l + {6m + f)yT'^] + 5R/2 



(16) 



X < 1 — exp 



5R/2 



l + (5«^ + /)V^'^ 

The number of nuclear fluorescence photons counted by the photodetector during 
time 9' between t'^ and tn+i is 

iV„ = iVe//-Pe.c(0-(l-e^^''), (17) 

where 

^Nr,-'-!^-^ (18) 

7 47r 

is the product of the number N^h of interrogated Thorium nuclei, the ratio '~fsn^/l of the 
rate of radiative decay of the isomer state to the total decay rate enhanced by possible 
electronic bridge processes, and the probability to detect the emitted photon. Here k is 
the quantum efficiency of the photodetector and Q is the effective solid angle covered 
by the detector. 

In turn, the excited state population at time t„+i is equal to 

at time t'^^^ it is equal to 

hR/2 



(19) 



Pexcifn 



+1) 



Pexc(Wl)eXp 
+ 



(7 + 

R 



i + (/-UVr'' 



7[i + (/-U7r'2] + 5i?/2 



(20) 
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and in turn 

iV„+i = iVe//-Pe(Ci)-(l-e-"'')- (21) 
Let us introduce the "evolution rate in the presence of laser field" 

and the functions 



a(A) = ^ ■ (1 - e-^' j ■ ^1 - ^ J ■ (1 - e-«(^)^) , (23) 

6(A) =e-^(^)^-^^'. (24) 
We then can rewrite expressions f fT7|) and fl2T]) as: 

iVn = a{f + 6m) + b{f + 6m)N^-i, (25) 

iV„+i = a{f - 6m) + b{f - 6m)Nn, (26) 

where A^„_i is the number of photon counts collected in the last measurement of the 
previous interrogation cycle, i.e. between t'^_^ and t„. 

We suppose that the generator offset / is small, i.e., / ^ T'. We then can expand 
the functions a and b near the working point 6m into a Taylor series: 

a{f ± 6m) =aoT aj, b{f ± 6m) = &o ± bj. (27) 
Using equations (125|) . ( 126|) and ( 127|) . we express the frequency offset / as 

. ^ Ar„+i -Nn- bojNn - iVn-l) . ^ 

^ 2ai-6i(iV„-iV„_i) ■ ^ ^ 

Now let us calculate the error 6f in determining the frequency offset /. Here we 
suppose that Nn+i and Nn are Poissonian random numbers with mean values described 
by fl25|l and fl26|l . i.e. the error of A'^; is a/A^. The average value Ni is 



N, ~ ao + boN, ^ (29) 

From f l25|) — f l29|) we express the error of the frequency offset as 

^ Vao{l-bo){l + bo + bl) ^^^^ 
a/2 [ai(l - bo) - biao] 
Substituting explicit expressions for oq, bo, ai, and bi into f l5U]) . we obtain 

iT' + 6l/TrVT^^ 



5f 



2R6mA/^N,fj{l-e~^'') 



^(1 - (1 - e-«^-^^')(l + e-«^-^^' + e-2Ge-27e') 
^(1 - e-Ge)(l - e-G^-^'^') + (!-§) ^e-Ge(l - e"^^') ' 

where G = G{6m)- 
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Expression fl3ip represents a fundamental lower limit on the error of the frequency 
offset for one interrogation cycle as depicted in Figure |H If the generator has good 
long-term stability, and the algorithm for frequency stabilization succeeds in the total 
correction of the measured offset (!28l) . then expression (13T]) is the frequency error 
corresponding to an absolute frequency instability of the solid-state nuclear clock after 
one interrogation cycle. 

The frequency error fl3T]) depends on the specific properties of the quantum 
discriminator {Nxh, 7, ^'), the detection system (kQ), on the intensity of laser radiation 
determining the excitation rate R, and on the regime of interrogation {6, 6', and 6m)- 
In practice the total interrogation time t = 2{9 + 9') will be imposed by the short-term 
stability of the interrogation laser (i.e. the duration over which frequency drifts and 
fluctuations are smaller than 5f). Also the laser intensity (and therefore the excitation 
rate R) will have a technical upper limit. On the other hand, the working point 5m and 
the relative irradiation time 



can be adjusted freely. It is therefore interesting to determine the optimal values 
of tR and 5m and study the dependence of the fractional frequency error 5f /u as a 
function of total interrogation time t and excitation rate R for optimized t^ and 5m- 
Because the general expression fl3T|l seems to be too complex for a direct interpretation, 
we will first focus on limiting cases of short and long interrogation times without 
restricting to explicit experimental parameters. In section [5] we will present numerically 
optimized results for the specific system of Thorium-doped Calcium fluoride crystals. A 
comparison between the analytic expressions and numerical optimizations can be found 
in Appendix B. 

4-2. Frequency error for short interrogation time 

Short interrogation time t = 2{9 + 9') refers to 76*' ^ 1 and G9 ^ 1, so the irradiation 
time and the detection time are shorter than the overall decay time and the inverse 
excitation rate respectively. In this case we can expand the exponents in (13T1) into a 
Taylor series and keep only the main terms. We obtain 



the absolute frequency error scales as with the total duration of the interrogation 
cycle. First we fix 5m and t and find the optimal value of t^. Therefore we should 
minimize the function 



5/ 



iX' + e/rP^ V3 {G - 7) {G9 + 7g"0 



(33) 



2R5m ^/^N:rf {0 + 0') v/yW 



GtR + 7(1 - tR) 

tnil - tR) 
producing the optimal value for tR. 




(34) 



^-7 
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The frequency error 6f for the optimized value of t^j is 

^"i?((5„,/r')v/5iv:;7' 

Taking into account f l2^ and minimizing ([35]) we can find optimal values for the working 
point 6m and the frequency error 6f. 

In the limit of weak field, i.e., for i? <^ 7, we find 

r' 

Sm = and (36) 

9r' 

^ ^t^l^^/RNlfi ^^^^ 
In the limit of strong field, i.e., if i? S> 7, we find 

5m = r', and (38) 



30 r' 

V = (39) 

4.3. Frequency error for long interrogation time 

We now consider the case where the interrogation time is long compared to the total 
decay time, i.e. when 7^ ^ 1. Here we can not neglect the exponents in (I3T!) because 
the optimal value either of 6 or of 6' can be comparable to I/7 or 1/G respectively. 
From the structure of fl3Tl) we expect a small optimal value of 6 because of the presence 
of y/l — e~^^ in the numenator. Therefore we suppose that 7^' ^ 1 and neglect the 
terms in e"'''^'. Under these asumptions the frequency error of the nuclear clock can be 
expressed as 



G(r + e/nV(i-e-)(i-«) 

2i?(5„^/5iV7;7(a(l - e"^) + (1 - a) a;e"^) 

where x = GO and a = '-//G. To find the optimal value for 6 we minimize (HOl) with 
respect to x which leads to the transcendental equation for x: 

e^-l = (41) 

a + 2(l-a)(l-x) ■ ^ ' 

This equation yields the optimal value of 9 = x/G for the case of long total interrogation 
time. We again consider the cases of weak and strong interrogation field. 

If the laser field is weak, i.e., if i? <^ 7, a is close to 1. Then 1 — a = 
5i?/(27(l + Si/T'^)) and dH]) transforms to 
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This equation contains a large term proportional to 7/-R on the right side of the 
denominator which has to be compensated by another large term 2xq — 1 to yield a 
solution xq of ( 14T|) . We obtain 



^o-^fl + Hl+l- (42) 



5R V r'2 

Therefore, for the case of weak laser field G9 ~ xq ^ 1 and we can neglect e"*^^ in ( HOj) . 
Optimization for 6m yields 

r' 

= and (43) 



We now consider the case where the laser field is strong, i.e. when i? ^ 7. We then 
can set a = in f HTj) . We obtain: 
xo 0.643798 

^=G=^^- ^^^^ 
In this case the optimal working point 6m and the frequency fluctuation 6f are: 

5m = r', and (46) 
Sf='-^. (47) 

The last expression determines the ultimate stability limit that can in principle be 
attained by the fluorescence spectroscopy method in a single interrogation cycle for a 
given quantum discriminator. It is valid for a large range of parameters as discussed 
in Appendix B. The optimized value for 6f ( H7|) is exclusively determined by the 
decoherence rate T' (which is estimated in Appendix A) and the number of interrogated 
nuclei, it does not depend on the exact values of the isomer transition energy or lifetime. 
To reach this performance level, the conditions i? <^ F' and R <^ Tggi have to be fulfilled, 
otherwise an accumulation of populations in quantum states not interacting with the 
laser field will decrease of Nefj over time. However, the very large ratio of r'/7 ~ 10^ 
gives a broad range of parameters where i? S> 7 and i? -C F', Tggi can be satisfied 
simultaneously. 

Note that in the limit of long interrogation time, the equilibrium value of the excited 
state population is attained in a short fraction of the interrogation cycle and the optimal 
excitation time takes only a small fraction of the interrogation cycle. In contrast, in 
the limit of short interrogation time, the equilibrium value can not be attained in one 
interrogation cycle and the optimal excitation time 6 should be equal to the photon 
collection time 6'. 
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In this section we present numerically optimized relative irradiation times t/j, ideal 
working points 6m, and hence ultimate frequency performances of the solid-state nuclear 
clock scheme for different excitation rates R and interrogation times t. We use the 
specific system parameters for ^^^Thorium ions doped into CaF2 crystals introduced 
in the previous sections: 7 = 7^?^-^ = 7.8 ■ 10~^ s~^, T' = 150 Hz ■ 27r = 942.5 s~^, 
u = 27r-1900THz = 1.19-10^^ s~^. We suppose that the total number of Thorium nuclei 
is Nth = 10^^, the solid angle is = An/lO and the quantum efficiency of photodetector 
is = 0.1 (which is typical for photomultipliers in our wavelength range). We obtain 
N^ff = 10^^ see ([18]). 

The optimized parameters 6m and tn are represented in Figure [5] for different 
excitation rates R. One can see that for all interrogation times, the optimal working 
points 6m he between T' and T' / \/2 which are the limiting results for _R <^ 7 f l5B| 
and > 7 dsn |46]) repectively. 

In general, the ideal working point is close to the inverse decoherence rate of the 
order kilohertz (see Appendix A). This means that the linewidth of the interrogation 
laser has to be significantly narrower than this, otherwise feedback stabilization to the 
nuclear transition will fail. Also one can see that the optimal value of t^ decreases 
with decreasing t and increasing R which means a shorter excitation than detection 
time. This is connected with the fact that the optimal illumination time 6 should not 
significantly exceed l/G{R,6m), i-e. the typical timescale to approach the equilibrium 
population of the excited state = |(1 — j/G{R, 6m)) according to f|T^ . 

Let us now turn to the nuclear clock error per one interrogation cycle obtained 
for the optimized parameters 6m and tji, expressed as fractional frequency error in 



a) 



10- 



^•7 

10-2 0.1 1 




10 



R = 0.037 

--R = O.OOI7 



r' 



r' 



10'' 



10'^ 



V2 



b) 



10- 




t-7 

0.1 1 



10 



R = O.OOI7 
R = 0.037 
i? = 7 
R = 3O7 
R = IOOO7 



10 10= 103 
i (s) 



10= 



Figure 5. Dependence of the optimal working point Sm (a) and of tji (b) on the total 
interrogation time t for different values of the excitation rate R. 




Figure 6. Fractional frequency error juj of the nuclear clock for different excitation 
rates R and different total interrogation times t. 

Figure [6l For short interrogation times, i.e., t < 0.1 7"^ the fractional frequency error 
scales as 5f /to oc t~^^'^ in accordance with (|37l) and ( !39l) . For long interrogation times, 
i.e. t > 107"^, the frequency error 6f/u becomes independent of the interrogation 
time in accordance with and fjTTj) . We note also that for strong pumping, i.e. 
R > 30 7, the clock performance is practically independent of the exact value of R, 
for long interrogation times this is already the case for i? > 7. For t > IO/7 and R 
more than a few 7, the fractional frequency error 6f/u per one interrogation cycle is 
essentialy independent of R and t and is equal to the ultimate value of about 1.8 ■ 10~^^ 
in correspondance with f H7|) . 

We now consider the parameter Sf/tu- \/t which describes the fractional clock 
instability. It should be optimized with respect to interrogation time t if one wants to 
perform a precise measurement of the frequency of the isomer transition (in the solid- 
state environment) by a long-term series of interrogations. It is also the relevant quantity 
describing the instability of the solid-state nuclear frequency standard when a narrow 
band laser is feedback-stabilized to the transition and compared to a primary frequency 
standard. If the series of measurements takes a time t = N ■ t, then the fractional 
instability ay^r) of the series is equal to the fractional frequency error averaged over the 
series. In shot noise limit we obtain 

(ry{T) = —y= = —j=. (48) 

Hence, to reduce Oyij) one should optimize (minimize) 5f\/t/uj = o-yy/r. We find 
that the optimal value for t is of the order of a few 7"^, only weakly depending on 
the excitation rate R, see Figure [7](a). The ultimate clock instability is given by the 
minimum value of Cyy/r at the optimal value of t. It saturates at ~ 1.63 x 10"^'^ s^/^ 
for high excitation rates (-R/7 ^ 10) and scales as for lower excitation rates 

{R/j ^ 0.1), see Figure [7](b). The resonance intensity Iq of the laser source required 
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Figure 7. Nuclear clock fractional instability Oyyfr = 5f/uj-^/tas& funtion of 
total interrogation time t for different excitation rates R (a) and the ultimate clock 
instability cr°P*y^ obtained for optimal value of t, as a function of R (b). 



to attain the "high excitation rate regime" -R/7 > 10 can be estimated using f lT^ . For 
the system parameters used in this section it yields Jq — 16mW/cm^. Other optimized 
parameters corresponding to i? = IO7 are: t ^ 8.88/7 — ^ ■ 10"^ s, t^ ~ 0.026 and 
5m ^ 0.991'. 

A few words on the vahdity of approximation (Q and the optimal value of the 
temperature: taking R = IO7 and using (1151] . we easily estimate the "reasonable" value 
27T ■ 0.53 Hz. The mixing rate Tggr should be large in comparison with 

27r ■ (5 — 100) Hz for room 
i (5 • 10-5 - 10-3) -T^IK]. 



of Veg as Veg ~ 



99 



this value. Our estimation (see section 12. 2p gives F 
temperature, and Tggi oc T^, then we can take Fgg'[Hz] 



Condition ^ gives T[K] > 



0.53 



or T > (20 - 100) K. If Q is not fulfilled, 



510-5-10-3 ' 

in a real multilevel system a significant fraction of nuclei will be pumped into the 
dark state, and the growing of the excitation rate R will reduce the performance. 
On the other hand, even at liquid nitrogen temperatures the contribution of the 
temperature-dependent 2nd order Doppler broadening term to the total decoherence 
rate F' becomes comparable to the temperature- independent magnetic dipole interaction 
term, and increasing the temperature is undesirable. We expect that the optimal sample 
temperature should be around the liquid nitrogen temperature. A more detailed analysis 
requires more precises data on the quadrupole relaxation and is beyond the scope of this 
article. Nevertheless we can claim that the optimal excitation rate R should be of the 
order of a few 7. 

In an experimental implementation, the performance of the solid-state nuclear clock 
will be limited by technical constrains such as the temperature stability of the sample, 
the frequency and intensity stability of the interrogation laser system, and noise of the 
detection system. If these limitations can be overcome, the solid-state approach can 
reach the 10"^^ fractional instability level within a few hours of averaging. It may hence 
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be competitive with the trapped-ion scheme [T7] concerning stabiUty for a (temperature 
stabihzed) hquid nitrogen temperature sample. Note that the fluorescence stabilization 
scheme developped here is opperated continuously without deadtime between cycles 
as no state initialization or re-trapping of particles is required. The global line shifts 
(several hundreds of MHz) due to crystal fields will obviously degrade the clock accuracy 
and require calibration with a primary standard. 

6. Conclusion 

We have analyzed the possibility of constructing a "solid-state nuclear clock" from 
^^^Thorium nuclei implanted into VUV transparent Calcium fluoride crystals. We expect 
static global line shifts to be dominated by the temperature-dependent electric monopole 
and electric quadrupole interactions and by 2""^ order Doppler shifts, necessitating active 
temperature stabilization of the sample. The electric quadrupole interaction leads to 
a splitting of the nuclear ground and isomer levels of the order of a few hundred MHz 
and is the main source of state mixing. As broadening effects we considered again the 
2nd Qj~(jgj~ Doppler shift, the electric quadrupole, and the magnetic dipole interaction. 
The magnetic dipole interaction does not depend on temperature and is hence expected 
to be the dominating source of decoherence (followed by the 2"^^ order Doppler effect) 
at liquid Nitrogen temperatures. For the CaF2 crystal we estimate a decoherence rate 
of 150 Hz. The extreme gap between this rate and the nuclear state lifetime renders 
the application of commonly used coherent interrogation schemes (i.e. Rabi, Ramsey) 
inappropriate. We therefore developed an alternative clock stabilization scheme based 
on fluorescence spectroscopy. A detailed analysis reveals a favorable scaling with 
the number of interogated nuclei, validating the solid-state approach. Neglecting all 
possible technical constrains (i.e. frequency and intensity stability of the interrogation 
laser), supposing that the sample contains 10^^ Thorium nuclei and that the detection 
probability is 1%, we find an ultimate fractional instability limit of 1.8-10"^^ for a single 
interrogation cycle taking about 10^ s. 
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Appendix A. Decoherence and mixing due to magnetic dipole interaction 

A Calcium fluoride crystal is formed from ^^F Fluorine and a natural mixture of Calcium 
isotopes (^°Ca, ^^Ca, ^^Ca, ^^Ca, and ^^Ca). The only Ca isotope with non-zero nuclear 
magnetic moment ^'^Ca has a small abundance of the order of 0.14 % whereas the nucleus 
of the only abundant ^^F isotope has a finite magnetic moment hf = 2.63 fiN and a 
nuclear spin s = 1/2. In the following we consider only the magnetic moment of the 
^^F nuclei in the CaF2 lattice to interact with the Thorium nuclear transition. The 
relaxation of the nuclear magnetic moment of Thorium (or Fluorine) nuclei is caused 
by the interaction with a random magnetic field "H at the position of the considered 
nucleus, created by the surrounding Fluorine nuclei. 

For the calculation of the relaxation and mixing rates we use the formalism of 
relaxation theory presented in [18]. Here we briefly describe this formalism and apply 
it first to the magnetic moment of Fluorine nuclei. This yields the characteristic rate 
of Fluorine spin-spin relaxation and therefore the characteristic correlation time of "H. 
Then we calculate the relaxation rate for the ultraviolet transitions between different 
quadrupole sublevels of the ground and isomer state of the Thorium nucleus. 



Relaxation method 



Here we follow a method developed in [18] for the calculation of the effect of a stationary 
random perturbation on the elements of the density matrix. The evolution of the density 
matrix is described by the Liouville-von Neumann equation 
dp i 



dt 



h 



H,p , where H = Hq + + Hi. 



(49) 



Here -B] = AB — BA is a commutator, = ^Ylal'^) '^o'i'^l Hamiltonian 
of the unperturbed system, H^' is the Hamiltonian of interaction with the coherent 
external field and Hi is the Hamiltonian of interaction with a random field responsible 
for the relaxation, in our case with the random magnetic field "H. In order to 
eliminate the strong contribution Hq from P9l) one applies the unitary transformation 
L ^ L' = Q'^^ot/hj^^-iHot/h the operators in fH^ (interaction representation). We 

obtain the equation for the density matrix: 
dp' 



dt 



I 

h 



Hi ,p 



I 

h 



dt'. 



The solution of this equation can be written as 

P'it)=p'{0)-'-J^ [{H^' + H'i),p'{t' 

Substituting (15T!) into the second term on the right-hand side of (!50!) yields: 
dp' 



(50) 



(51) 



dt 



I 

n 
1 



H^'{t),p'{t) 



% 

h L 



H',{t),p'{0) 



H'.it), iHt'{t') + H',{t')),p'{t 



dt'. 



(52) 
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We assume that the characteristic time scale of the correlations of fluctuation of the 
random field H[{t) is much smaller than all other characteristic time scales, including 
the characteristic time scale of variation of the density matrix [§|. Then one can choose 
the origin in such a way that the density matrix would not change noticeable during this 
time interval (one may substitute p'{t') by p'{t)) but all the correlations of the random 
field will disappear (one may take — oo as a lower limit of integration). Also suppose 
that all the matrix elements of H[{t) have zero mean values (otherwise we would add 
them to Hi')- Performing the averaging over the ensemble (indicated by overline), we 
obtain 

^^'=4[^f'w,p'(t)]-^ 



dt 



H^'{t),p'{t) 



H[it), H[{t')rp\t) 



dt'. 



(53) 



The first term in fl53|) corresponds to the evolution of the density matrix under the action 
of a coherent field. The second term describes the relaxation. If the perturbations are 
stationary, it can be written as 

1 



5'(t) 



f9 



H[{t), H[{t-e),p'{t) 



d9. 



Expanding the commutators in (1541) . one obtains: 



Ua'P')t], 



(54) 



(55) 



where ojj 



R„ 



jal3a'l3'{^ap) + 3 P' a' Pa{^ P' a') 



1 

7 7 ^ 



and 



Hiap{t)Hip,o,'it - 0) e 



de 



(56) 



(57) 



(note that in the right-hand side of the last expressions, the matrix elements of Hi are 
taken in the "initial" Schrodinger representation). Expression (l55l) for the relaxation 
operator R contains oscillating terms that can be neglected. As a result we keep only 
the terms where Uap = ^a'p'- Except for some special cases, this condition is fulfilled 
only if a = a' and f3 = /3' (mixing of populations), or ii a = (3 and a' = (3' (relaxation 
of coherences). 

Using (j54l) — (j57l) one can calculate the relaxation and population mixing rates. 



§ Although this condition is not entirely fulfilled for the relaxation of Fluorine spins in the random 
magnetic field % created by neighboring spins, quite good correspondancc with experimental results 
indicates the applicability of this theory. For the Thorium nucleus interacting with this condition 
is fulfilled much better due to the much lower magnetic moment compared to Fluorine. 
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Spin relaxation of a Fluorine nucleus 

Let us consider a single Fluorine nucleus. The spin of this nucleus can take one of two 
projections m = ±1/2 with respect to the arbitrary quantization axis z. Let us denote 
|1) = |m = —1/2) and |2) = |m = 1/2). The magnetic moment operator fl for the spin 
of this nucleus can be expressed as /2 = sfip/s. Therefore we can express the operator 
Hi = —fl ■ % as 

i^i = /iF(^.(|l)(l|-|2)(2|) 

-?/,.(|l)(2| + |2)(l|)-zH,(|l)(2|-|2)(l|)). (58) 

In this work we use the "diffusion" model for the correllations of the random 
magnetic field: 

nt)) = 0, (ii\t)) = nl mt)nj{t - 9)) = 5.,^e-l^l/*^ (59) 

The average square of the random magnetic field "Hq and the correlation time tc will be 
determined later. 

According to fl56|) the relaxation rate -R1212 of the coherence pi2 is 

^1212 = ^ bn22(0) + j221l(0) - j212l(0) 

- J2222 (0) - Jim (0) - J2121 (0)] . (60) 
From (l57l) and (159|) one finds 

^1212 = ^ [jll22(0) + J221l(0) 

- J2222(0) - Jllll(O) - 2j212l(0)] . (61) 

where 

Jllll(O) = J2222(0) = -Jll22(0) = -J221l(0) = ^^t^ 

j2i2i(0)=2^t, (62) 

which yields the coherence relaxation rate 

8 fipT-Lltc . s 

Ruu - -3-^- (63) 

Similarily one finds the population mixing rate as 

-Riiii = ~ -R1122 

J2121 (0) + jl212 

3 ' ^ ^ 

Now we can estimate the correlation time tc which is determined by the population 
relaxation rate: 
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Figure 8. Notation for Tliorium nucleus spin states used in this paragraph. 

To find "Hq we recall that in CaF2 the Fluorine nuclei form a simple qubic lattice with a 
distance of a/2 = 2.73 A between two neighbors. Spins of different Fluorine nuclei are 
not correlated, therefore the average square "Hq of the random magnetic field H is the 
sum of squares of magnetic fields produced by different Fluorine nuclei: 

nl = Y,K- (66) 

a 

T-La is the magnetic field created by a-th Fluorine nucleus whose radius- vector is with 
respect to the nucleus of interest. Using jl^ = 2sa^F, and Ha = {?>r{jlj^a) — f^a''^a) / 
we find 

-8 

a 'a 



Therefore 



= ^ E ['^ + "^ + "T •8.4 (68) 

l,m,n=—oo " 

(here the sum is taken over l,m,n excluding / = m = n = 0), or "Hq ^ 21.5 which 
leads to a relaxation time tc — 15 /is according to (I65p . Note that this value is quite 
close to experimentally observed spin-spin relaxation times ~ 20 /is |49j . 

Spin relaxation of a Thorium nucleus 

We now consider an individual Thorium nucleus interacting with random magnetic 
fields created by the surrounding Fluorine nuclei. The quantization axis is fixed by the 
direction of the main axis of the electric field gradient. Let us denote the different states 
of the Thorium nucleus as \i), where i is a number from 1 to 10 corresponding the spin 
state shown in Figure [81 The average square T-Lxh magnetic field at the position 

of the Thorium nucleus is equal to 



E f E (*-^) 



6/i 



F 



^6 



■19^48.5Gl (69) 

\ z / 

l,m,n=~oo '- k=l,m.,n 

The relaxation rates Fjj = —Rijij for optical coherences pij due to interaction with 
the random magnetic field created by surrounding Fluorine spins can be calculated in 




Figure 9. (a): dependence of Sf ■ t^l"^ for short interrogation time t on the excitation 
rate R obtained by the optimization of (PT|) (curve 1), using approximate formula ([55)) 
(Unc 2), and (|37p (hue 3). (b): dependence of 5/ for long t on R obtained by the 
optimization of ([3T|) (curve 1), using approximate formula (|47)) (line 2), and (|44|) (line 
3). 



the same manner as for a single Fluorine spin above. These rates (for allowed magnetic 
dipole transitions) are equal to: 

'H%htc 19 



r4,8 

r2,8 
Taj 



r3,9 = 

r4,9 = 
r4,io 



7 



r 



6,10 



2 f 



2 19 2' 



fx 



9 

25' 
19 



9 



- 27r ■ 142 Hz, 

- 2n ■ 150 Hz, 
27r-84 Hz, 

- 27r- 158 Hz, 

- 27r ■ 108 Hz, 



3h? 



5' ' 25' 
[/ig - 2/ie/ig + /i^] ~ 27r ■ 251 Hz, 



(70) 



Here we used the values of the Thorium nuclear magnetic moment in the ground state 
Hg = OAbfiN and in the isomer state fie = —0.076fiN respectively [50]. The decoherence 
rates Fa § = F4 9 have been used in the evaluations of the nuclear clock performance in 
section |5l 



Appendix B. Validity of analytic expressions for the frequency fluctuation 

Here we compare the numerically optimized results for 6f obtained in section [5] to 
the analytic expressions obtained in section H] for the limiting cases of short and long 
interrogation times. For short t, 5f oc t"^/^, therefore we express as a function 
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of R, and compare theses results with the results of the analytical formulas (1371) and 
fl39|) . see Figure [9](a). For long t we express 6f as a function of R and compare this 
results with formulas (l44l) and (147|) . 

We conclude that the approximate analytic formulas (137|) . (l39l) . (144|) . and (1471) 
well describe the behavior of the optimized frequency error 6f of the nuclear clock for 
different limiting cases of short and long interrogation time and for weak and strong laser 
fields. Taking into account the quite smooth dependence of Sf on the interrogation time 
(see Figure Ej), we conclude that one can use expressions ( 1371) . ( l39l) . ( 14^ . and ( 147|) as a 
lower-bound estimation for the frequency fluctuation 6f of the fluorescence spectroscopy 
interrogation scheme. 
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